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Abstract. This paper surveys our recent research on quantum infor- 
mation processing by nuclear magnetic resonance (NMR) spectroscopy. 
We begin with a geometric introduction to the NMR of an ensemble of 
indistinguishable spins, and then show how this geometric interpretation 
is contained within an algebra of multispin product operators. This al- 
gebra is used throughout the rest of the paper to demonstrate that it 
provides a facile framework within which to study quantum information 
processing more generally. The implementation of quantum algorithms 
by NMR depends upon the availability of special kinds of mixed states, 
called pseudo-pure states, and we consider a number of different meth- 
ods for preparing these states, along with analyses of how they scale 
with the number of spins. The quantum-mechanical nature of processes 
involving such macroscopic pseudo-pure states also is a matter of debate, 
and in order to discuss this issue in concrete terms we present the results 
of NMR experiments which constitute a macroscopic analogue Hardy's 
paradox. Finally, a detailed product operator description is given of re- 
cent NMR experiments which demonstrate a three-bit quantum error 
correcting code, using field gradients to implement a precisely-known 
decoherence model. 
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available as vol. 3385 from the International Society for Optical Engineering, 1000 
20th St., Bellingham, WA 98225, USA. 



1 Introduction 



It has recently proven possible to perform simple quantum computations by 
liquid-state NMR spectroscopy [12, 13, 14, 15, 16, 22, 32, 42]. This unprece- 
dented level of coherent control promises to be quite useful not only in demon- 
strating the validity of many of the basic ideas behind quantum information 
processing, but more importantly, in providing researchers in the field with new 
physical insights and concrete problems to study. This is particularly true since 
the ensemble nature of the systems used for NMR computing differs substan- 
tially from the systems previously considered as candidate quantum computers. 
The use of ensembles provides tremendous redundancy, which makes computa- 
tion with them relatively resistant to errors. It also has the potential to provide 
access to a limited form of massive classical parallelism [14] , which could for ex- 
ample be used to speed up searches with Grover's algorithm by a constant but 
very large factor [9, 24, 36]. The barriers that have been encountered in extend- 
ing NMR computing to nontrivial problems further raise interesting questions 
regarding the relations between microscopic and macroscopic order, and between 
the quantum and classical worlds [23, 43]. 

NMR computing is also contributing to quantum information processing 
through the assimilation of theoretical and experimental NMR techniques. These 
techniques have been developed over half a century of intensive research, and 
grown so advanced that a recent book on the subject is entitled "Spin Chore- 
ography" [21]. It is noteworthy that, due to the scope of its applications, NMR 
is now more often studied in chemistry and even biochemistry than it is in 
physics and engineering, where it was initially developed. This has had the ef- 
fect that a large portion of these techniques have been discovered empirically 
and put into the form of intuitive graphical or algebraic rules, rather than devel- 
oped mathematically from well-defined principles. Thus the interest which NMR 
computing is attracting from the quantum information processing side likewise 
has the potential to benefit the field of NMR spectroscopy, particularly through 
the application of algorithmic, information theoretic and algebraic techniques. 

Finally, NMR has the potential to contribute in significant ways to the de- 
velopment of its own mathematics, in the same ways that computers have con- 
tributed to the development of recursive function theory, number theory, com- 
binatorics and many other areas of mathematics. By performing experiments 
which can be interpreted as computations in homomorphic images on the alge- 
bras that are naturally associated with NMR spectroscopy, it may be possible to 
obtain insights into, or even "proofs" of, algebraic properties that would other- 
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wise be inaccessible. To give some idea of its potential computational power, we 
point out that the spin dynamics of a crystal of calcium fluoride one millimeter 
on a side, which can be highly polarized, superbly controlled and measured in 
microscopic detail by NMR techniques [53, 57], is described by an exponential 
map in an algebra on about 4 1 " 11 physically distinct dimensions. 

This paper will survey our recent research on quantum information process- 
ing by liquid-state NMR spectroscopy, including some new experiments which 
serve to clarify the underlying principles. We begin with a geometric interpreta- 
tion of the quantum mechanical states and operators of an ensemble of identical 
spin 1/2 particles, both pure and mixed, which provides considerable insight 
into NMR. The corresponding geometric algebra is then extended to the product 
operator formalism, which is widely used in analyzing NMR experiments, and 
which constitutes a facile framework within which to study quantum informa- 
tion processing more generally [48] . We proceed to use this formalism to give an 
overview of the basic ideas behind ensemble quantum computing by liquid-state 
NMR spectroscopy, with emphasis on "pseudo-pure" state preparation and scal- 
ing. Next, we consider one way in which quantum correlations can appear to 
be present even in weakly polarized spin ensembles, and illustrate this with the 
results of NMR experiments which constitute a macroscopic analogue of Hardy's 
paradox [26] . Finally, the utility of NMR and its associated product operator for- 
malism as a means of studying decoherence will be demonstrated by an analysis 
of our recent experiments with a three-bit quantum error correcting code [15] . 4 



2 The geometry of spin states and operators 

NMR spectroscopy is based on the fact that the nuclei in many kinds of atoms 
are endowed with an intrinsic angular momentum, the properties of which are 
determined by an integer or half-integer quantum number S > 0, called the 
nuclear spin. For the purposes of quantum information processing by NMR, it 
will suffice to restrict ourselves to spin S = 1/2. In this case, measurement of 
the component of the angular momentum along a given axis in space always 
yields one of two possible values: ±h/2 (where Ti is Planck's constant h over 2tt). 

4 The reader is assumed throughout to be familiar with the basic notions of quantum 
information processing, as presented in e.g. Refs. [43, 44, 51, 56]. Excellent detailed 
expositions of NMR spectroscopy are also available, see e.g. Refs. [19, 21, 40, 47]. A 
more introductory account of our work on ensemble quantum computing by NMR 
spectroscopy, directed primarily towards physicists, may be found in Ref. [16]. 
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According to the principles of quantum mechanics, the quantum state of the 
"spin" (nucleus) after such a measurement may be completely characterized by 
one of two orthonormal vectors in a two-dimensional Hilbert (complex vector) 
space H, with Hermitian (sesquilinear) inner product (-|-). A rotation of this 
axis in physical space induces a transformation in H by an element of the special 
unitary group SU(2), which is the two-fold universal covering group of the three- 
dimensional Euclidean rotation group SO (3), and the elements of TC are called 
spinors to emphasize this fact. The Lie algebra basis (I x ,I y ,I z ) of SU(2) (or 
S0(3)) corresponding to infinitesmal rotations about three orthogonal axes in 
space satisfies the commutation relations 

[J x ,J y ] = ll z , [I Z ,I X ] = lly , [Iy,Iz] = tlx , (1) 

and the eigenvalues ±1/2 of these three Hermitian (self-adjoint) operators cor- 
respond to the possible outcomes of measurements of the angular momentum 
(in units of K) along the three axes. 5 

The Hilbert space representation of the kinematics of an isolated spin, how- 
ever, is not sufficient to describe the joint state of the macroscopic collections 
of spins which are the subject of NMR spectroscopy. A mixed state (as opposed 
to the pure state of an isolated spin) is a random ensemble of spins not all in 
the same pure state. This "ensemble" could be a thought-construction which 
describes our state-of-knowledge of a single spin (as used in J. W. Gibbs' for- 
mulation of statistical thermodynamics), or it could be a very large physical 
collection of spins, as in an NMR sample tube. In either probability is 

assigned to every possible spinor, which can be interpreted as its frequency of 
occurrence in the ensemble (but see Ref. [30] for a Baysian point-of-view) . The 
Heisenberg uncertainty principle limits what can be known about the ensemble 
to the ensemble-average expectation values of the quantum mechanical observ- 
ables. This information, in turn, can be encoded into a single operator on H, 
called the density operator. 

To define this operator mathematically, we first recall the canonical algebra 
isomorphism between the endomorphisms End(7i) and the tensor product Tt ® 
7i* of 7i with its dual space H*. Denoting the dual of a vector under 
the Hermitian inner product of Ti by (ip |, the composition product in End(7i) 

5 Detailed explanations of these basic features of the quantum mechanics of spin may 
be found in modern textbooks. We would particularly recommend Sakurai [45], for 
an introduction to the underlying physics, or the monograph by Biedenharn and 
Louck [5], for a complete mathematical development. 
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corresponds to a product on H <£> Ti* which is given on the factorizable tensors 

by 

= {<p'\#){\<p)®(#\), (2) 

and extended to all tensors by linearity. Following common practice, we shall usu- 
ally drop the tensor product sign "(g)" and write this dyadic product as \ip)(i)\. 
The restriction of this product to the diagonal, | ip ) ( ip \ , linearly spans the (real) 
subspacc of all Hermitian operators in End(H), and the action of SU(2) on these 
products is its usual action on such operators, 

IVXVI ~ u\i>){il>\u , (3) 

where U = U~ denotes the Hermitian conjugate (adjoint) of U £ SU(2). 

Restricting ourselves to an ensemble involving a finite set of states {| ipk}} 
for ease of presentation, the density operator may now be defined as [7] 

P = \^}(lP\ = J2 k Pk\^k)(^k\ , (4) 

where the pu > are the probabilities of the various states in the ensemble 
(J2kPk = !)• Because (<p\p\tp) = J^kPk \{v\^k)\ 2 > for any spinor \ip), 
the density operator is necessarily positive semi-definite. Letting "tr" be the 
contraction operation on Ti ® Ti* (or trace on End(W)), letting A £ End(7i) 
be any Hermitian operator, and using the invariance of the trace under cyclic 
permutations, we find that 

tr(Ap) = ^2 k Pk tr ( A | Vfc > < ^fe | ) = J2 k Pk ( ^ ' A ' ^ k ) ' ( 5 ) 

This proves our claim that all ensemble-average expectation values can be ob- 
tained from p. Note in particular that tr(p) = 1. 

Before showing how this applies to NMR spectroscopy, we wish to introduce 
an important geometric interpretation of the spin 1/2 density operator, and in- 
deed of the entire operator algebra. As operators, the angular momentum com- 
ponents transform under SU(2) by conjugation, i.e. I w UI w U (w £ {x, y, z}). 
Thus U and —U induce the same transformation, so that conjugation consti- 
tutes a transitive group action of S0(3) on the real linear space ( I x ,I y , I z ). It 
follows that this space is naturally regarded as a three-dimensional Euclidean 
vector space. Note further that ( 1, J x , J y , I z ) equals the four-dimensional space 
of Hermitian operators on Ti, where 1 is the identity on Tt which we will hence- 
forth identify with the scalar identity 1. This shows that any density operator 
can be uniquely expanded as the sum of a scalar and a vector: 

p = ±tr(p) + tr(/ x p)2/ x + tr(/ y p)2/ y + tr(/ z p)2/ z = ±(1 + P) (6) 
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We call P the polarization vector, since its length P = \\P\\ < 1 (the polariza- 
tion) is a measure of the overall degree of alignment of the spins in the ensemble 
along P. The positive semi-definiteness of p implies P < 1, and if P = 1, the 
density operator describes a (ensemble of spins in the same) pure state up to an 
overall phase factor. In this latter case the density operator can be written as 
p = | ip)(il)\ for some spinor and hence is idempotent (equal to its square). 
This vectorial interpretation of two-state quantum systems became widely known 
through the work of Feynman, Vernon and Helwarth [20] , although it is inherent 
in the phenomenological equations for NMR first proposed by F. Bloch [6] (see 
below). 

To extend this geometric interpretation to the entire algebra generated by 
( 1, J x , J y , J z ), we regard the composition product of angular momentum opera- 
tors as an associative bilinear product of vectors. We shall call this the geometric 
vector product. Since the eigenvalues of the (S = 1/2) angular momentum op- 
erators are ±1/2, the eigenvalues of their squares are both 1/4, from which it 
follows that (2J X ) 2 = (2J y ) 2 = (27 z ) 2 = 1. In accord with the isotropy of space, 
moreover, 

(UI w Uf = U{I w fU = 1/4 (7) 

for all U € SU(2) and w <G {x, y,z}, which together with the bilinearity of the 
product implies that the square of any vector is equal (relative to the orthonor- 
mal basis (27 x , 2J y , 2J Z )) to its length squared. Via the law of cosines, we can 
now show that the symmetric part of the geometric product of any two vectors 
is their usual Euclidean inner product: 

A B =\(U\\ 2 + \\B\\ 2 -\\A-Br) 

= \ (A 2 + B 2 — (A - B) 2 ) = \ (AB + BA) 

The commutation relations in Eq. (1), on the other hand, show that the anti- 
symmetric part is equal (up to a factor of —i) to the usual vector cross product: 

Ax B = -%[A,B] = -\{AB — BA) = -i(AAB) (9) 

We call the antisymmetric part A A B = i{A x B) the outer product of A 
and B, and note that it is geometrically distinct from vectors because inversion 
in the origin does not change it. Such things have been called "axial vectors", 
although wc prefer the older and more descriptive term bivector. On writing 
the geometric product as the sum of its symmetric and antisymmetric parts, 
AB = A- B + A/\B,we see that perpendicular pairs of vectors anticommute. 
It follows that the three basis bivectors i2I x , i2I y and i2I z also anticommute. 
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These square to —1 rather than 1, however, and thus can be identified with 
the usual quaternion units [1]. Finally, the unit pseudo-scalar 8I x I y I z likewise 
squares to — 1, which together with the fact that it commutes with the basis 
vectors and hence everything in the algebra enables it to be identified with the 
unit imaginary i itself [25]. 

This algebra is often called the Clifford algebra of a three-dimensional Eu- 
clidean vector space, although we shall use the term geometric algebra here 
(which W. K. Clifford himself used). Such an algebra is canonically associated 
with any metric vector space, and provides a natural algebraic encoding of the ge- 
ometric properties of that space. The fact that the three-dimensional Euclidean 
version can be defined starting from the well-known properties of the spin 1/2 
angular momentum operators indicates that a large part of quantum mechanics 
is really just an unfamiliar (but extremely elegant and facile [27, 29]) means 
of doing Euclidean geometry. Geometric algebra has more recently been exten- 
sively advocated and used to demystify quantum physics by a number of groups 
[4, 17, 28, 38]. Of particular interest are recent proposals to use the geometric 
algebra of a direct sum of copies of Minkowski space-time to obtain a relativistic 
multiparticle theory, from which all the nonrelativistic theory used in this paper 
falls out naturally as a quotient subalgebra [18]. 

We are now ready to describe the simplest possible NMR experiment. The 
time-independent Schrodinger equation is 

%h\j>) = H\tp) , (10) 

where the Hamiltonian H is the generator of motion and the "dot" denotes the 
time derivative. This implies that the density operator evolves according to the 
Liouville-von Neumann equation: 

ifip = ihY,kPk (\^k)(ipk\ + \^k)(^k\ S j 

= EkPk(H\^k)(i>k\-\A)(A\H) = [H,p] 

The dominant Hamiltonian in NMR is the Zeeman interaction of the magnetic 
dipoles of the spins (which is parallel to their angular momentum vectors) with 
a constant applied magnetic field B . This Zeeman Hamiltonian is given by 
Hz = —7}jhB 0} where 7 is a proportionality constant called the gyromagnetic 
ratio, which together with the above gives the Bloch equation [6]: 

P = p = -%^[p,B ] = jPxB Q (12) 

The solution to this equation is p(t) = Up(0) U with U = cxp(— liifz), which 
is a time-dependent rotation of the polarization vector about the magnetic field 
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with a constant angular velocity ujo = 7^11 Bo II- This "classical" picture is an 
example of Ehrenfest's theorem, and is analogous to the precession of a gyroscope 
in a gravitational field. 

Throughout this paper we adopt the universal convention that the magnetic 
field is along the z-axis: B = B 2I Z . The component of the net precessing 
magnetic moment of the spins in the transverse xy-plane generates a complex- 
valued radio-frequency electrical signal proportional to tr((J x + il y )p(t)) = 
2J X • P(t) + i2I y ■ P(t), whose Fourier transform is an NMR spectrum con- 
taining a peak at the precession frequency of each distinct kind of spin present 
in the sample. This has the important consequence that in NMR we measure the 
expectation values of the observables directly, which is due in turn to the fact 
that we are measuring the sum of the responses of the spins over the ensemble. 
These measurements yield negligible information on the quantum state of the 
individual spins in the ensemble and hence are nonperturbing, in that they do 
not appreciably change the state of the ensemble as a whole. Such weak mea- 
surements contrast starkly with the strong measurements usually considered in 
quantum mechanics, where determining the component of a spin along an axis 
yields one of two possible values and "collapses" it into the corresponding ba- 
sis state, so that only one classical bit of information can be obtained [43, 45]. 
A discussion of the computational implications of weak measurements may be 
found in Ref. [14]. 

The natural (minimum energy) orientation of the spins' dipolcs in a magnetic 
field is parallel to the field, and thus to obtain a precessing magnetic dipole it is 
necessary to rotate the polarization vector P away from the field axis 21 z . This 
is done by applying an additional, rotating magnetic field B\ of magnitude B\ in 
the xy-plane perpendicular to the static field B , which gives the time-dependent 
Hamiltonian 

H = Hz + H R f = -jh(B I z + B 1 (cos(ujt)I x + sm(ujt)Iy)) . (13) 

The effect of such a rotating field is most easily determined by transforming 
everything into a frame which rotates along with it, in which the Hamiltonian 
becomes time-independent: 

p' = e'^'pe 1 ^ 1 ' , H' = e-^titHe 1 "* 1 * = Hz + ^hBrl^ (14) 

Then the Bloch equation itself is transformed as follows: 

P' = -iujI z P' + e -*»ti*p e ™ti* + PW Z 

(15) 

= P' x {H 1 - ojI z ) 
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Thus if ijj equals the natural precession frequency of the spins uio = jhB , the 
Zecman Hamiltonian H' z = Hz = uqI z cancels out. In this frame, the spins 
turn about a (rotating) axis perpendicular to B' x at a rate uj\ = jfrBi, so that if 
the polarization vector starts out along z, it is in the xy-plane where it produces 
the maximum signal after a time t = 7r/(2wi). Henceforth, all our coordinate 
frames will be rotating at the transmitter frequency unless otherwise mentioned. 

3 The product operator formalism 

Thus far we have restricted our presentation to ensembles consisting of indistin- 
guishable nuclear spins. The power of NMR spectroscopy as a means of chemical 
analysis, however, depends on the fact that the different nuclei in a molecule gen- 
erally have distinct electronic environments, which affect the applied magnetic 
field at each nucleus. As a result, they precess at slightly different frequencies 
and give rise to resolvable "peaks" in the resulting spectrum. This is also one 
of the reasons why NMR provides a facile approach to quantum information 
processing, since it permits each chemical equivalence class of spins in the en- 
semble to be treated as a separate "qubit". In this section we will describe an 
extension of the density operator to multispin systems, using a basis which is a 
direct generalization of the "scalar + vector" basis given above for a single spin. 
We then illustrate this so-called product operator formalism [8, 19, 48, 50, 54] 
by describing how quantum information processing can be done on an ensemble 
of multispin molecules, using the internal Hamiltonian of liquid-state NMR. For 
the sake of simplicity we shall assume throughout that the ensemble is in a pure 
state, i.e. that the joint state of the spins in every molecule is identical. The next 
section is devoted to the complications involved in extending this approach to 
the highly mixed states which are available in practice. 

As usual in quantum information processing [51, 56], we choose a computa- 
tional basis (| ), | 1 )) for the Hilbcrt space H of each spin that corresponds to 
the eigenvectors of its J z operator, i.e. to alignment of the spin with (up) and 
against (down) a magnetic field Bo along the z axis. Relative to this basis, a 
superposition co|0) + Ci| 1) (c ,ci ^ complex with |c | 2 + |ci| 2 = 1) is any 
state with transverse (xy) components. The Hilbert space needed to describe the 
kinematics of a system consisting of N distinguishable spins (not an ensemble) 
is the AT- fold tensor product of their constituent Hilbert spaces [43, 45]. The 
induced basis in this (2 JV )-dimensional space is 

Ik 1 )® \k 2 }®---® \k n ) ee \k 1 k 2 ...k n ) EE |fc), (16) 
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where K n € {0,1} (n = 1, . . . , N) is the binary expansion of the integer k e 
{0, . . . , 2^ — 1}. Because of the canonical isomorphism 

End(W) ® End(H) ss End(H ® W) (17) 

together with our previous isomorphism End(7Y) w H H* , this implies that 
the density operators for an ensemble of iV-spin molecules are all contained in 
the iV-fold tensor product space 

{H <8> • • • ® H) <8> (W* <8> • • • <8> W*) « (ft ® W) <8> • • • ® (W <g> W) • (18) 
It follows that a basis for the algebra of TV-spin operators is 

\k)(£\ ^|« 1 « 2 ...^)(A 1 A 2 ...A W | 

= (| K 1 )(A 1 |)®(I« 2 )(A 2 |)®---®(|« A, )(A W |), 

where k", A™ e {0, 1} (n = 1, . . . , AT) are binary expansions of the integers k, I e 
{0, . . . , 2 W — 1}. This basis, however, does not consist of Hermitian operators, 
and although the dyadic products 1 1\) ) ( ip \ (\i>) € H ® • • • ® H) do span the real 
subspace of all Hermitian operators, the restriction of the basis in Eq. (19) to 
the diagonal does not. 

An algebra basis which has the advantage of also being a linear basis for the 
subspace of Hermitian operators is known as the product operator basis. It is 
induced by the one-spin basis (1, J x , J y , J z ), and consists simply of the tensor 
products of the angular momentum operators of the individual spins. In the case 
of two spins, this basis has sixteen elements: 

10 1 1(g) J x 1(g) J y l(gJ z 

J X (g 1 J X (g J X J X (g Jy J X (g J Z ^ 
Jy (g 1 Jy (g J X Jy (g Jy Jy <g> J Z 
J z (g 1 J z (g J x J z (g Jy J z (g J z 

As before, a notation which eliminates the need for repetitive "<g" symbols is 
preferred. This is obtained by using superscripts for the spin indices in the op- 
erators 

J™ = l(g---(gl(gJ„,(gl(g---(gl (21) 

(I w in the n-th place, n = 1, . . . , N, w € {x, y, z}), and noting that by the mixed 
product formula between the operator composition and tensor products: 

I™I n v = l(g---(gJ tt (gl(g---(gl(gJ, ) (g---(gl = J^J™ (22) 
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(J„ in the m-th place, I v in the n-th, m,n — 1, . . . , N with m < n, and u, v e 
{x, y, z}). In the following, we will again identify the identity operator 1 <g> ■ • ■ (g> 1 
with the scalar identity 1. We will also be using the operator norm ||-f^,!| 2 = 
((J^,) 2 ) = (I^) 2 = 1/4 obtained from the scalar part, rather than the more 
usual Frobenius norm ji = tr((J^,) 2 ) = 2 N ~ 2 on End(W), because the former 
is independent of N. The normalization of our basis to || J™|| = 1/2 rather than 
1 will be seen to have both advantages and disadvantages, but the convention is 
well-established in NMR. 

Just as with an ensemble consisting of a single type of spin, a pure state 
may be characterized by the idempotence of its density operator: p 2 = p. The 
scalar part of the density operator is (p) = 2 _Jv tr(p) = 2~ N , and if we write 
an arbitrary density operator p = 1 1\> ) ( | in diagonal form as 

/ 2 W -1 \ 

P = U {E k=0 Pk\k)(k\jU (23) 

(0 < Pk < 1) J2kPk = 1) f° r some U e 511(2^), we see that the idempotence of 
p is equivalent to ( p 2 ) = 2~ N , i.e. pe — 1 for some I G {0, . . . , 2 N — 1}. This 
shows that the density operator of a pure state is in fact a primitive idempotent. 
Without loss of generality we may take I = 0, so that p = U\ 00 . . . )( 00 . . . \U. 
If we expand | )( | in the product operator basis, we obtain E + = |(1 + 2J Z ), 
and similarly for 1 1 )( 1 1 = E_ = |(1 - 2J Z ). Thus we can also write the density 
operator of a pure state as 

p = U (E\ E\ ■ ■ ■ if , (24) 

where the superscript on the idempotent E + is the spin index as before. More 
generally, the set of all density operators consists of the closed convex cone of 
positive semi-definite operators in the Hermitian subspace of End(W<g) ■ • • ® H), 
and the density operators of pure states are the extreme rays of this cone. 

We now consider the form of the Hamiltonian which is operative among 
the spins of an ensemble of molecules in the liquid state (with which we are 
exclusively concerned in this paper) , again using the product operator formalism. 
First, there is the Zecman Hamiltonian previously given for a single spin, i.e. 

H Z ee -WqJ* ^ I? (25) 

with ujq = 7i7 n (l — a n )Bo, where 7™ the gyromagnetic ratio of the n-th spin 
and < a n < 1 is the chemical shift due to the (usually small) influence of 
the electronic environment of the spins on their precession frequencies. This 
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Hamiltonian is easily seen to be diagonal in the computational basis \k) (k = 
0, . . . , 2 N - 1), with eigenvalues (±w* ± ■ ■ ■ ± uj$ r )/2. 

Second, there is an exchange interaction known as the J or scalar coupling, 
which is proportional to the inner product of the spins' polarization vectors, 
namely 

H } = V 2ir J mn (J™I™ + I™ I" + 7™I?) , (26) 

where J m ™ is the coupling strength in Hertz. This interaction is mediated by 
the electrons in the chemical bonds between atoms, and is usually negligible for 
atoms separated by more than three bonds. Standard perturbation theory [45] 
shows that the eigenvalues of the total Hamiltonian H = Hz + H j are given to 
first order by the diagonal elements of 

H' = H z + H'j = -Hz + 2ttV J m ™/™7™ , (27) 

whereas the eigenvectors are given to first order by: 

w - l'>^-E m ..,^E w ^r^w^r (28) 

The numerator of each term in the summations is nonzero only if n p = X p for 
p 7^ m, n and K m = (1 — A"), in which case it is irJ mn , while the denominators 
of the corresponding terms are — ujq . It follows that the eigenvectors are 
negligibly perturbed so long as the frequency differences are much larger than 
the scalar couplings, i.e. |w™ — » 7r| J mn \. We shall be making this weak 
coupling approximation throughout. 

Another, potentially quite large term in the molecular spin Hamiltonian is a 
through-space interaction between the spins' magnetic dipoles. Because of the 
rapid motions of the molecules in a liquid, however, these interactions are aver- 
aged to zero much more quickly than they can have any net effect. The effective 
absence of this interaction nevertheless has the important consequence that the 
spins in different molecules do not interact, and hence cannot be correlated with 
one another. 6 As a result, the density operator of the entire sample g (which 
describes an abstract Gibbs ensemble obtained by tracing over the spins' envi- 
ronment) can be factorized into a product of density operators for the individual 
molecules, i.e. 

Q = Q 1 ■ ■ ■ Q M = P 1 ® • • • ® P M , (29) 

6 More precisely, the spins do not interact to an excellent, but first-order, approxi- 
mation; second-order effects do exist and are a source of spin-spin relaxation (aka 
decoherence) . 
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where g m = 1 ® • • • ® p m (8) • • • ® 1. In a pure liquid (or if we are looking at just 
one component of a solution), all the molecules are equivalent so that all these 
density operators are the same. It follows that we can work with the partial 
trace over all but any one of the molecules, which is called the reduced density 
operator p = p 1 = ■ ■ ■ = p M . Since this operates on a space of dimension 2 N 
where N is now the number of spins in a single molecule, rather than 2 MN where 
M ~ 10 20 is the number of molecules in the sample, this is a very considerable 
simplification. It also means that in liquid-state NMR we are working with a 
physical ensemble (the sample), rather than a purely abstract Gibbs ensemble. 

Finally, there is the interaction of the spins with a transverse RF (radio- 
frequency) field, which we described in the last section. Whenever weak coupling 
is valid, we can apply this field in a single "pulse" , tuned to the precession 
frequency of the fc-th spin (say), which is short enough that we may neglect the 
evolution of the spins due to scalar coupling while it lasts. This effects a unitary 
transformation of the form 

e-^ = l- t (f)2/«-I(f) 2 + Hf) 3 2/x + 
= cos (|) -tsin(|) 2l\ , 

which corresponds to a right-hand rotation of the fc-th spin by an angle 9 about 
the x axis in the rotating frame. Using a pulse with a broad frequency range, it is 
also possible (in fact easier) to apply such a rotation to all the spins in parallel. 

We will now indicate how RF pulses, in combination with the innate Hamil- 
tonian of the spins, enable us to implement standard quantum logic gates in a 
manner similar to that considered by computer scientists studying universality 
in quantum computation [3]. The simplest such gate is the NOT operation on 
the e.g. first spin, which simply rotates it by 7r; combining the above formula 
with the basic geometric algebra relations I X I Z = I Z I X and (2J X ) 2 = 1, we 
obtain 

e ^il E i + ^il = { -2ill)El(2zll) = |(1 + 8Ji 

^(l-S/ 1 / 1 / 1 ) = ±(1-2/1) = E i_ 

The c-NOT (controlled-NOT) gate, on the other hand, is a n rotation of e.g. 
the first spin conditional on the polarization of a second. Using the relation 
E\E 2 -^ — 0, we can easily show that 

{-2il\E 2 _ + E 2 + ){E\E 2 ± ){2illE 2 _ +E 2 + ) = E\ e E 2 ± (32) 

(e e {±})- The phase factor i multiplying l\ complicates the action of the c- 
NOT on a superposition, but can be eliminated by a phase shift conditional on 
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(30) 



the second spin. Using E 2 __ + E 2 + = 1, this phase-corrected c-NOT gate is given 

2I 1 X E 2 _ +E 2 + = {-%E 2 _ + E 2 + )(2iI x E 2 _ + E 2 + ) 

= (l + (e" l § - 1)£ 2 _) (l + (e' wJ * - l)E 2 _^j (33) 

_ e -tf E 2 _ e nrllE 2 _ _ e -i§(l~2ll)E 2 _ 

and hence the idempotents J5± also give us an algebraic description of the c- 
NOT gate, in addition to the density operators of pure states. It is well-known 
that single spin rotations, together with the c-NOT, arc sufficient to implement 
any quantum logic gate [3]. 

The c-NOT can be implemented in NMR by combining single spin rotations 
with the conditional rotations induced by (weak) scalar coupling 2i:J 12 l\l 2 
[16, 22, 31]. Recalling that in the discrete S0(3) subgroup of rotations by n/2, 

e-'i'yVi 1 * 1 = e'i^e'i 1 * 1 , (34) 
we can expand the above propagator as follows: 



S^ 2 = e^i ( 1 - 2I l) E '- 



■,f I] e -znElE 2 _ g*§ I) 



(35) 



Since the overall phase of the transformation has no effect on the density oper- 
ator, it follows that the c-NOT gate S 1 ' 2 can be implemented by an NMR pulse 
sequence, wherein each pulse and delay corresponds to the indicated "effective" 
Hamiltonian in temporal order: 

[-f 7 y ] [nil I\\ -+ [-§ 1\\ -+ [— § (l\ + II)] 

exp (f {l\ + I 2 )) cxp (f I\) cxp {-*l\l 2 ) cxp (fly 1 ) 

In practice, the effective Hamiltonian [irl\l 2 ] is obtained by applying a 7r-pulse 
to both spins in the middle and at the end of a 1/(2 J 12 ) evolution period, to 
"refocus" their Zeeman evolution [16]. The [— f-fy] and [— Hamiltonians are 
implemented by RF pulses as above, while the [— § (J z + I 2 )] transformation is 
most easily implemented by letting one spin evolve while applying a 7r-pulse to 
the other, then vice versa, and finally realigning the transmitter's phase with 
the spins'. 
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The "readout" procedure needed to determine the result of an NMR compu- 
tation differs somewhat that usually considered in quantum computing [51, 56]. 
The most important difference is of course the fact that in conventional NMR 
one can only make weak (nonperturbing) measurements of the observables, as 
previously described. As likwise described above, these observables are the x and 
y components 7™ and 7™ of the dipolar magnetization due to each spin in a ro- 
tating frame defined by the transmitter frequency. The products of the angular 
momentum components of different spins (e.g. 7 x 7 y ), however, do not produce a 
net magnetic dipole and hence cannot be detected directly. Thus we are limited 
to one-spin observables, as is usually assumed in quantum computation. The un- 
observable degrees of freedom may also be characterized in the basis | k ) ( I | as 
having a coherence order | ( fc | 7 Z | fc ) — ( ^ | J z | ^ ) | 7^ 1, where 7 Z = l\ + • — Y I z 
is the total angular momentum along z [16, 19, 21, 40, 47]. 

According to the usual phase conventions of NMR, the Fourier transform of 
the x-magnetization of e.g. the first spin, l\, yields an absorptive peak shape, 
while 7 y produces a dispersive shape, both centered on its precession frequency 
ujq. If the first spin is coupled to e.g. the second, its signal is modulated by 
cos(7rJ 12 i) yielding a spectrum containing two peaks separated by the coupling 
constant J 12 [16]. An effective exception to the unobservability of the products 
are those of the form (or 7 y 7 2 ), which (when J 12 7^ 0) evolve under scalar 
coupling into one-spin terms. Using the facts that 47* I 2 and 47* l\ anticommute 
while (47 z 7 2 ) = 1, we can show this as follows: 

e - l t7rJ 12 4/ z 1 / z 2 = e -*t2TrJ 12 l\ll (4/1/2) e4t2 7rJ 12 / z 1 / z 2 

= (cos(tt J 12 t) - iAl\ll sin(^ J 12 t)) (47i7 z 2 ) (37) 
= cos(tt J 12 t) + sin(7r J 12 t) 27 y 

Because the signal is now sinusoidally modulated by the coupling, for a single 
pair of coupled spins this results in a pair of antiphase peaks with opposite signs, 
as opposed to the inphase peaks described for 7* and 7 y above. These antiphase 
peaks may likewise be absorptive (7*7 2 ) or dispersive (7 y 7 2 ), respectively. Figure 
1 shows examples of all these possibilities for a pair of two coupled spins. 

More generally, if the n-th spin is coupled to M others, its signal is split into 
2 M peaks at frequencies of (u^/tt ± J mi " ± ■ ■ ■ ± J mM ") /2, one for each combi- 
nation of "up" and "down" states for the M spins to which it is coupled. If the 
transverse magnetization is due to a ir/2 rotation of a spin polarized along z as 
before, then the heights of these peaks are proportional to the probability differ- 
ences between pairs of states | n mi . . . n n . . . n mM ) <-> | K mi ... (1 — k") . . . K mM ) 
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Inphase Inphase Antiphase Antiphase 

Absorptive Dispersive Absorptive Dispersive 




Fig. 1. Plots of NMR spectra for a weakly coupled two-spin molecule (amplitude versus 
frequency). The left-hand plot is for the spin state which gives a pair of inphase 

absorptive peaks for spin 1 (left) and a pair of inphase dispersive peaks for spin 2 
(right). The right-hand plot is for the spin state + illy, which gives a pair of 

antiphase absorptive peaks for spin 1 (left) and a pair of antiphase dispersive peaks for 
spin 2 (right). By fitting the peak shapes in such spectra after various -k/2 rotations 
of the individual spins, one obtains sufficient information to uniquely reconstruct the 
complete density operator. 



separated by flips of that spin. To show this, we restrict ourselves to two spins 
for ease of presentation, and consider a general diagonal density operator of the 
form 



Pzz =Po|00)<00|+pi|01)<01|+p2|10)<10|+p3|ll)<ll| 

= \ + s(Po + Pi - P2 - Pa)l\ + ^(Po-Pi+P2- Ps)ll (38) 
+ (Po -Pi -Pi + p 3 )l\l\ , 



where pk denotes the probability that a molecule is in the state \k). Rotating 
this to 



= i + \ijPo + Pi - P2 - Pa)ll + ^{po-Pi+P2- p 3 )ll ( 39 ) 

+ (PO -Pi -P2+P3,)l\ll 
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and computing the signal in the Zccman frame yields 

tr (e-^^'^pxz e^ J "'^ (ij + i/J)) 
= |((Po + Pi -P2 - p 3 ) cos(7rJ 12 t) ^ 

+ (Po - Pi - Vi + P3) * sin(7r J 12 t)) 
= ie- jl2 *(po -P2) + \e~™ J "\vx ~Ps) , 

thus showing that the peaks at ujq ± it J 12 have amplitudes proportional to the 
probability differences as claimed. 

In closing, we mention that although vector interpretations of single quantum 
inphase (J x ) and antiphase (ill 2 ) states are available (and widely used in NMR 
[21]), no satisfactory geometric interpretation of general product states is known. 
The development of an intuitive model for the geometry determined by the action 
of 511(2^) on the product operators thus stands as an open problem in the field. 
There are two reasons why the problem is nontrivial. The first is the well-known 
the existence of correlated states, whose density operators cannot be factorized; 
in the case of a pure state, these states are also called entangled [43]. We shall 
consider such states further in Section 5. The second, much less widely recognized 
reason is that there is but one imaginary unit for all the spins, so that in the 
tensor product of their geometric algebras the unit pseudo-scalars 87™/™/™ must 
be identified by taking an appropriate quotient [48]. This is a form of implicit 
correlation which is always present even in otherwise factorizable states. Further 
discussion of this issue may be found in Refs. [17, 18]. 



4 Pseudo-pure state preparation and scaling 

Liquid state NMR must be done at temperatures far above the differences be- 
tween the spin Hamiltonian's energy levels (eigenvalues). The ensemble's spin 
state thus represents a compromise between the constant force of the applied 
magnetic field and the forces of the random fields induced by the thermal mo- 
tions of spins in other molecules. Thus pure states are not available, so that the 
underlying ensemble is not uniquely determined by its density operator. This 
would seem to make NMR useless as a means of performing deterministic com- 
putations, but in fact a class of mixed states has been found for which a state 
vector is (up to an overall phase) canonically associated with the density op- 
erator. This section is devoted to describing the properties and preparation of 



17 



such pseudo-pure states, with emphasis on the computationally important issue 
of how they scale with the number of spins N. 

According to the principles of quantum statistical mechanics [7], the den- 
sity operator p eq for an ensemble of TV-spin molecules at thermal equilibrium 
is given by the Boltzman operator determined by their common Hamiltonian, 
exp(— H/k B T), divided by the corresponding partition function Z eq = 
tr(exp(— H/ksT)) (where ks is Boltzman's constant). The Hamiltonian H is 

N tN 
i z ■ 



well-approximated by its dominant Zeeman term Hz 
Given the gyromagnetic ratios of nuclear spins and the strongest available mag- 
netic fields, we have ||i!fz||/(fcBT') <~ 10~ 5 at the temperatures needed for liquid- 
state NMR, so that a linear approximation is quite accurate: 



-"eq 



1 - H z /k B T 
tr(l - H z /k B T) 



1 - H z /k B T 
2 N 



(41) 



In homonuclcar (i.e. single spin isotope) systems, one can assume that ujq = 
hB {l — a n )-f n rts HBo'j is constant for all n. Since the amplitude of an NMR 
signal is also proportional to imprecisely known factors determined by the spec- 
trometer setup, Wq/^bT) is usually set to unity when analyzing a homonuclear 
experiment (or to the ratios of each 7™ with min m (7 m ) otherwise). The parti- 
tion function 2~ N is likewise constant for any given system, but because of our 
interest in scaling we shall always include it explicitly in this section. 

It is important to observe that, because the angular momentum components 
observed by NMR have no scalar part (i.e. are traceless), the scalar part (identity 
component) of the density operator 2~ N does not contribute to the signal. It also 
does not evolve under unitary transformations, and hence NMR spectroscopists 
usually forget about it altogether — even though it comprises the vast majority of 
the norm of the density operator. In these terms, the equilibrium density operator 
of a two-spin system, and its matrix representation in the usual computational 
basis, is 

" -1(11+11) = 



±Diag(l, 0,0,-1) 



00>( 00 I - 1 11 >< 11 1 

"1 0" 





-1 



(42) 



where the "hat" on p eq signifies its traceless part. 

In contrast, the density operator of two spins in their pseudo-pure ground 
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(assuming 7 > as usual) state is 

Poo - ±W\ + II + Kill) = ±§ (E\E\ - \) (43) 
= ±|(|00)<00|-|) <-> ±^Diag(3, -1,-1,-1). 

The overall sign depends on whether we have a population excess or deficit in 
the ground state; for consistency, we shall generally assume the former. Observe 
that a unitary transformation of the density operator induces a transformation 
of the corresponding state vector just as it does for true pure states, since 

Up oa U = !(([/|00»(E/|00>r-i) . (44) 

Similarly because the NMR observables A = J™, J™ are traceless, the ensemble- 
average expectation value relative to a pseudo-pure density operator yields the 
ordinary expectation value versus the corresponding state vector: 

tr(Ap 00 ) = ±(tr(A|00)(00|)-±tr(A)) = §<00|A|00) (45) 

The general form of a pseudo-pure density operator is 

Pi, = ^hW^\~ 2 ~ N ) > ( 46 ) 

where | ip ) is a normalized TV-spin state vector, and the prefactor has been cho- 
sen so as to keep the maximum eigenvalue ||p^||2 equal to that of the iV-spin 
equilibrium density operator. 

Even though we have defined them to have the same maximum eigenvalue, 
for iV > 1 the remaining eigenvalues of p eq and differ, and hence there is no 
unitary transformation taking one to the other. There are nevertheless a number 
of nonunitary processes by which one can prepare pseudo-pure states. The most 
direct is to generate a spatially varying distribution of states across the sample, 
such that the ensemble average is pseudo-pure. This can be done by using a 
field gradient along the z-axis to create a position-dependent phase shift whose 
average is zero, thereby in effect setting the transverse (xy) components of the 
density operator to zero. 7 For example, it is readily shown that the sequence 

[f (J* + II)] - [nlll 2 z ] - + 7 y 2 )] (47) 

applied to the two-spin equilibrium state p eq yields 

2-1 (V3 (E\El \ Illl) llEl Bit) , (48) 

7 In the homonuclear case, the zero-quantum coherences are not rapidly dephased by 
a z- gradient, so a slightly more complicated procedure is necessary. 
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which is reduced by a z-gradient to (3/32)5 {E\E\ - 1/4). Further RF and 
gradient pulse sequences which convert the equilibrium state of two and three 
spin systems to pseudo-pure states may be found in Ref. [16]. 

An alterative proposed by E. Knill et al. [34] is to "time- average" the results 
of several separate experiments. In the simple case of a two spin system, the 
average of the three states 

P123 = + ~ ±Diag(l, 0,0,-1) 

P23i = + | Diag(l, 0,-1,0) (49) 

p 312 = ±(2/ z / z 2 + / z 2 )^ ±Diag(l, -1,0,0) 

is the pseudo-pure state 

|(Pl23 + P231+P312) - ^(2J Z 1 + 2I Z 2 +4I Z 1 J Z 2 ) q 

<-» ^Diag(3, -1,-1,-1). 

More generally, one can obtain the same results that one would get on a pseudo- 
pure state by averaging the results of the experiments over all 2^ — 1 cyclic 
permutations of the nonground state populations of the equilibrium density op- 
erator. Although this naive approach is not efficient, Knill et al. have shown that 
one can average over smaller groups in time 0(N 3 ) with much the same effect. 

A fundamentally different approach, first proposed by Stoll, Vega & Vaughan 
[52] and subsequently adapted to NMR computing by Gershenfeld & Chuang 
[22], involves working with subpopulations of molecules distinguished by the 
states of additional ancilla spins. Gershenfeld and Chuang [12] have given an 
example of a two-spin conditional pseudo-pure state (as we call it), which is 
obtained by row/column permutation of the diagonal equilibrium density matrix 
Diag(3, 1, 1, —1, 1, —1, —1, — 3)/16 of a three-spin system including one ancilla, 
namely 

^Diag(3, -1,-1, -1,-3, 1,1,1) 

~ i(2/ z (2/ 2 + 2/ z 3 + 4/ 2 / z 3 )) (51) 

= \{E\-E\_){ElE%-\). 

The last form makes it clear that in the subpopulation with the first spin "up" , 
which is labeled by E + , and in the subpopulation with it "down", which is 
labeled by E_, spins 2 and 3 are in the pseudo-pure state E 2 + E^ + — 1/4. Since 
the spectrum of spins 2 and 3 is antiphase with respect to the ancilla spin 1, 
one can select the subpopulations just by keeping only either positive or negative 
peaks. Although the situation is considerably more complicated with more spins, 
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Gershcnfeld and Chuang have shown that conditional pure states can be obtained 
(with some loss of signal) using as few as 0(\og(N)) ancillac. 

An alternative to conditional pure states, which we call relative pseudo-pure 
states, can be obtained via the partial trace operation (in NMR, decoupling 
[19, 21, 40, 47]), rather than peak selection as above. For example, a two-spin 
relative pseudo-pure state is given by the partial trace over the ancilla spins 1 
& 2 in 

£ Diag(4, 2, 2, 0, 2, 0, -2, 0, 0, -2, 0, 2, 0, -2, -2, -4) 

<-» ^ (E\E 2 + {E% + E\) + E\E 2 _(E%E% - E 3 _E 4 + ) (52) 

+ E l _E 2 + (EtEt - E 3 + Ei) - E l _E 2 _{E s _ + E 4 _)) , 

which is again a permutation of the diagonal elements of p eq . This can be seen 
by adding up the 4x4 blocks of the matrix, obtaining Diag(6, —2, —2, — 2)/32. 
Alternatively, since the partial trace in the product operator formalism corre- 
sponds to simply eliminating those terms depending on cither spins 1 or 2 and 
multiplying the remaining terms by 4 [48], we need only add up the multipliers 
of E\E 2 + , . . . , E\E 2 _, which yields 

^((1 + E\ - E 3 _)(l + E\ Et) 1) 
= \ (E\E% - i) ~ ^ Diag(6, -2, -2, -2) . 

We now consider briefly how the SNR (signal-to-noise ratio) of these meth- 
ods of creating pseudo-pure states scales with the number of spins N. It has 
been argued that since the equilibrium population of the ground state falls off 
exponentially with the number of spins, and all these methods are aimed in some 
fashion at isolating the signal from the ground state population, the SNR of all 
these methods must likewise decline exponentially with N [55]. Although this 
argument carries considerable weight, we shall see that the number and variety 
of the available methods renders the actual situation rather more complex. The 
standard to which the signal strength must be compared is that of a single spin 
in its equilibrium state, namely 

P eq = \1\ = !(|0)<0|-|1)<1|). (54) 

The maximum eigenvalue ||p eq ||2 = 1/4 is what we will use as the standard 
signal strength for spins of like gyromagnetic ratio (as assumed throughout). 
We shall therefore calculate the SNR of a pseudo-pure state by transforming it 
to the corresponding ground state |0---0)(0---0| — 2~ N (if need be), taking 
the partial trace over all but one of the spins, and multiplying the maximum 
eigenvalue of the result by 4. 
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Fig. 2. Negative base ten logarithm of the polarization P as a function of the loga- 
rithm of the ratio of the energy level spacing to k B T for a one (solid), four (dash), 
sixteen (dash-dot) and sixty-four (dot) spin pseudo-pure state obtained by cyclic 
averaging. For protons in a standard 500 MHz spectrometer at room temperature, 
E/(k B T) « P ~ 10~ 5 at equilibrium. 

The maximum eigenvalue of the partial trace over all but one of the spins in 
a pseudo-pure state obtained by cyclic averaging, as in Eq. (50), is easily seen 
to be N/(4(2 N — 1)), which decays almost exponentially with the number of 
spins N. There is an additional factor of y/2 N — 1 which comes from averaging 
over 2^ — 1 experiments, and gives a net SNR of N/ (4y/2 N — 1) for the average. 
The exponential time requirements of cyclic averaging will nonetheless force one 
to average over smaller groups, with consequently smaller improvements in the 
SNR. In any case, the SNR declines superpolynomially with N. Figure 2 shows 
how the signal strength changes as a function of the ratio of the energy level 
spacing to k&T, relative to the signal in a perfectly polarized sample, when the 
pseudo-pure state is obtained by cyclic averaging, for varying numbers of spins. 

Because of the many possible variations on the ideas and the difficulty of 
analyzing all of them, it is not practical to present simple formulae for the SNR 
of the other methods of preparing pseudo-pure states. Further complexity is 
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added to the situation by the ability to combine the various methods above. A 
number of such combinations arc given in Knill et al. [34], along with bounds 
on the SNR for each. In our laboratory we are developing a new method, again 
based on field gradients, which enables the sample to be divided into discrete 
volumes and separate unitary transformations to be applied to each. In principle, 
this permits multiple experiments to be performed, and their results added, in a 
single experiment, thereby performing an average over multiple experiments in 
constant time. This new method could also be used in a variety of combinations 
with existing methods. 

It is nevertheless encouraging to observe that the SNR of the two-spin condi- 
tional and relative pseudo-pure states given in Eqs. (51) and (52) is 1/2 in both 
cases; this is exactly the decline in the ground state population of a two-spin 
system compared to a one-spin. In Eq. (51), we attain this "theoretical limit" be- 
cause the expansion of the density operator consists of a single term conditioned 
on the state of a single ancilla; it is not possible to do as well with more spins. In 
Eq. (52), however, it is because such permutations arc able to concentrate polar- 
ization in a subset of the spins. We have found this makes it possible to derive 
a two-spin pseudo-pure state from a six-spin equilibrium state with no loss of 
SNR, whereas a simplistic ground-state population argument implies we should 
lose at least 1/2. This may be seen by adding up the rows in the rearrangement 
of p eq shown in Eq. (55) below, which corresponds taking the traces of the four 
16 x 16 blocks along the diagonal, and yields Diag(48, —16, —16, —16). 

Diag(6, 4, 4, 4, 4, 4, 4, 2, 2, 2, 2, 2, 2, 2, 2, 2, (55) 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -2, -2, -2, -2, -4, -4, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -2, -2, -2, -2, -4, -4, 

2, 2, 2, 2, 2, 2,-2,-2,-2,-2,-2,-2,-2,-4,-4,-6 ) 

The partial trace over one of the two remaining spins then gives Diag(32, —32), 
which when divided by 128 (twice the partition function) yields \I T as claimed. 

A general algorithm has recently been given by Schulman & Vazirani [46] 
whereby one can "distill" an M-spin relative pure state from an ensemble of 
molecules each containing N spins. Starting from a uniform polarization of P, 
this algorithm yields M perfectly polarized spins providing M/N ~ 0(P 2 ), a 
result anticipated by earlier work in NMR which showed that the polarization of a 
single spin can be enhanced by at most a factor 0(y/~N) [49] . Unfortunately, given 
that P w 10~ 5 for protons at equilibrium in a standard 500 MHz spectrometer, 
a molecule with of order 10 10 spins would be needed to prepare a perfectly 
polarized state on a single spin — which is in a pseudo-pure state at equilibrium! 
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The importance of Schulman and Vazirani's algorithm thus lies in the fact that 
it shows that there is sufficient order in a typical NMR sample of 10 20 spins at 
room temperatures to make it at least theoretically possible to perform quantum 
computations on of order 10 10 spins. 

One might hope that a more tractable algorithm, in terms of the absolute re- 
sources required, could be found by requiring only that it produce a pseudo-pure 
state with bounded SNR from the high-temperature equilibrium state. Since in 
the high-temperature approximation the largest element of the density matrix 
decays exponentially with the number of spins, it is clear that any such an 
algorithm must go beyond that approximation. Even so, given that Schulman 
and Vazirani's algorithm is currently far beyond our ability to implement physi- 
cally, it seems unlikely that a practical breakthrough will be obtained by purely 
algorithmic means. Fortunately, physical methods of "refrigerating" the spins 
are available, for example optical pumping [41]. These are presently confined to 
very simple systems, but such a source of polarization could in principle be used 
in conjunction with polarization transfer techniques to produce (pseudo-)pure 
states on large numbers of spins. Even at the low polarizations we can con- 
veniently access, however, NMR has proved itself to be a powerful means of 
exploring quantum dynamics in Hilbert spaces of substantial size. To illustrate 
this, we will now present the results of NMR experiments which constitute a 
macroscopic analogue of a quantum mechanical test for quantum correlations 
that arc inconsistent with the existence of "hidden variables" . 



5 Macroscopic analogues of quantum correlations 

Given the success of the purely classical Bloch equations (and their multispin 
extensions) in describing liquid-state NMR phenomena [19, 21, 40, 47], it is 
perhaps surprising that experiments can be performed whose mathematical de- 
scription, at least, is formally identical to that of experiments which are believed 
to demonstrate uniquely "quantum" phenomena. For example, Seth Lloyd has 
recently proposed that the nonclassical correlations in (Mcrmin's version of) 
the GHZ state can be validated using NMR [37]. His approach involves using a 
fourth "observer" spin to perform a nondcmolition measurement on the three 
spins in a GHZ state (or a pseudo-pure analogue thereof). Here we shall describe 
experiments which demonstrate another, rather different way in which we can 
"emulate" quantum phenomena with liquid-state NMR. In reading this account, 
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it should be kept in mind that although pseudo-pure states do provide a faith- 
ful representation of the transformations of pure states within the highly mixed 
states that are available in liquid-state NMR, their physical interpretation dif- 
fers significantly from that of true pure states. Hence, as discussed in greater 
detail at the end of this section (and also in Ref. [11]), our results should not 
be taken to resolve any foundational issues in quantum mechanics [43]. They 
demonstrate, nonetheless, a degree of coherent control sufficient to enable such 
issues to be addressed, if these same transformations and measurements were 
applied to a true pure state. 

The approach taken here was inspired by an educational paper published 
a few years ago, in which T. F. Jordan has shown that the contradictions with 
hidden variables implied by violations of Bell's inequalities as well as by the GHZ 
and Hardy's paradox can be derived entirely by consideration of the expectation 
values of product operators, rather than by observations on single spins [33] . This 
shows that, in principle, it is not necessary to use nondemolition measurements 
with an observer spin in order to perform experiments which demonstrate these 
contradictions by NMR; it can be done directly from observations on ensembles 
of the spins of interest, providing at least they are in a (pseudo-)pure state. 
In a companion paper to Jordan's, N. D. Mermin points out that in real-life 
experiments it is nevertheless not possible to perform the measurements, either 
of single spins nor (by implication) of expectation values, with sufficient precision 
to establish the "perfect" (total) correlations on which "EPR" arguments against 
the existence of hidden variables are based [39]. In that same paper, however, 
Mermin shows that Hardy's paradox is a special case of the Clauser-Hornc form 
of Bell's inequality. This enables Hardy's paradox [10, 26] to be extended to 
an open set in the Hilbert space of only two spins, to which sufficiently precise 
experimental data can confine us. 

In the following, we present the results of NMR experiments which implement 
the specific example of Hardy's paradox presented by Mermin in an Appendix 
to his paper [39]. Let us map the "red" and "green" eigenstates 1 1G) and 1 1 R ) 
of Mermin's measurement 1 to the spin states | ) and 1 1 ) , respectively. It will 
be clearer here to relabel this measurement as "A" , and to use | «g ) = I ) 
and | aR ) = 1 1 ) as synonyms for its eigenbasis. Correspondingly, we will relabel 
Mermin's measurement 2 as "B", and denote its the eigenbasis by 

|&> = 7f|0>--/f|l) and |/3 R ) ee v /||0) + v /f|l). (56) 

Then the state which Mermin has shown leads to a near-maximum violation of 
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Bell's inequality while also providing an example of Hardy's paradox is 

|V> = i I 00 > + ./§ | 01 ) + ./§ 1 10 > . (57) 



To translate this into the context of NMR, we first note that the observable 
whose expectation value is the probability that measurement A yields the state 
| ctG ) is given by A = E + = ^(1 + 2J Z ) (we drop the usual spin index because 
the measurements A & B are assumed the same for both spins). Similarly the 
observable which gives the probability that measurement B yields | /3q ) is 

B ee | / 3 g) ( / 3 gH ||0)(0| + ||1)(1|- v /A(|0)(1| + |1)(0|) 

2 ^ 5 z V 25 

In addition, the density operator (including the identity) of Mermin's state is 



# EE \^)^\ = \ + \(l\+ll)-\l\ll 

+ yfj (11(1 + 211) + (1 + 2/1)/* 2 ) (59) 

+l(/^ x 2 +/;/ y 2 ) . 

The state | 01 ) is obviously related to | 00 ) by a rotation of spin 2, while | 00 ) 
can likewise be rotated to | 10), but without affecting |01), by a conditional 
rotation of spin 1. Wc shall denote these by 

P{(j>) ee e - '* 1 ? and Q(6) = e^ 67 ^ , (60) 

respectively. They act consecutively on the ground state to yield 

(00\P(4>)Q{6) = [cos(0/2)cos(<£/2), sin(6»/2) cos(0/2), sin(0/2), 0] , (61) 

which is easily verified to equal (ip \ = [1/2, y/3/8, ^/3/8, 0] when 

cj) = 2arctan(y375) and 6 = 2 arctan(y / 372) . (62) 



Using the product operator techniques presented in section 3, these transforma- 
tions are readily implemented by NMR pulse sequences. 

The next thing to notice is that if we take expectation values with the usual 
idempotents E\E 2 + , . . . , E X _E 2 _, we get 



I = 4(# E l + E 2 + ) ee A(&A 1 A 2 ) 

| =A{^E 1 + E 2 _)= 4 («P r A 1 (l - A 2 )) 

| = 4(^E 1 _E 2 + ) = 4 (^(1 — A 1 ) A 2 ) 

p2\ — a I iT r { -i /»2 



(63) 

= \i^E x _E 2 _) ee 4<«P(1 - A x )(\ - A 2 )) . 
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These correspond to the diagonal of the density matrix in the usual I z basis, 

diag(!F) = [|,|,f,0] (A on 1, A on 2), (64) 

which contains the probabilities of the four possible outcomes of performing 
measurement A on both spins (as shown). 

The product operator form of B immediately makes clear that measurement 
B is just a measurement of the magnetization of the spin along an axis inclined 
at an angle of ( = arctan(\/24) = 7r — 9 to the z-axis in the xz-plane. Letting 
R = exp(— iC,I y ), it follows that the probability that measurement B on spin 1 
yields "G" (i.e. | j3 G )) is 

4(«fB 1 ) = A{^R 1 A 1 R 1 ) = 4( J R 1 «P'i? 1 A 1 ) . (65) 

with a similar expression for spin 2. More generally, the probabilities of the 
outcomes of the other combinations of measurements are given by the diagonals 
of the transformed density matrices: 

diag (R 2 \PR 2 ) = [0, §, ^, ±] (A on 1, B on 2) 

diag (f? 1 ^ 1 ) = [0,^,1,^] (B on 1, A on 2) (66) 

diag (R'R^tfR 1 ) = [jjfe, |7 , i|] ( B on 1, B on 2) 

For compactness, let us denote these probabilities by tf^, where i, j € {A, B} are 
the measurements and k, I € {G, R} are the corresponding outcomes, e.g. 'Pq^ — 
4(*'A 1 (1-B 2 )). 

We may translate Mermin's proof [39] that these probabilities are incompat- 
ible with hidden variables associated with the individual spins into this context 
as follows: First, since &qq = \Pqq = 0, in any molecule wherein one of the spins 
is parallel to the z-axis the other must be antiparallel to the axis of measure- 
ment B and vice versa. Hence, since \Pqq is nonzero, in some molecules (9%, to 
be precise) both spins must be antiparallel to the z-axis. But this contradicts 
the fact that = 0. More generally, Mermin has shown that 

f G B G B < fG A G B + ^R A R A + ^G A (67) 

is an example of the Clauser-Horne form of Bell's inequality [43]. Hence to dis- 
prove the existence of such one-particle hidden variables it would be sufficient 
to determine these probabilities to ±2% or so. 

At this point we encounter a significant complication, which is that the 
"strong" (von Neumann) measurements assumed in their analyses by Jordan 
and Mermin cannot be implemented by NMR; we can only perform "weak" 
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Fig. 3. The pairs of 13 C-labeled chloroform (CHCi'a) spectra (carbon left, proton right) 
obtained by performing the four combinations of measurements A C A H , A C B H , B C A H and 
B C B H on the pseudo-pure form of Mermin's state & . The spectra have been normalized 
by the height of the peak of the corresponding spin in the pseudo-pure ground state, 
and the horizontal axis is in kHz. The transitions of the peaks, from left to right, are 
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Probabilities of Outcomes G & R for the Measurements 
A & B (Carbon, Proton) Demonstrating Hardy's Paradox, 
as Derived from the Chloroform NMR Spectra in Fig. 3 



Measurements 


(G,G) 


(G,R) 


(R,G) 


(R,R) 


Residuals 


(A C ,A H ) 


0.253 


0.380 


0.366 


0.001 


0.008 


(A C ,B H ) 


0.029 


0.609 


0.217 


0.145 


0.018 


(B C ,A H ) 


-0.002 


0.230 


0.614 


0.159 


0.005 


(B C ,B H ) 


0.097 


0.125 


0.156 


0.622 


0.021 
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(nonperturbing) measurements of the probability differences between states con- 
nected by single spin flips [16]. This is done by applying a magnetic field gradient 
along the z-axis, which (as previously described) dephases any transverse com- 
ponents in the density operator. Thereafter, a pair of "soft" n/2 readout pulses, 
each tuned to the frequency of just one of the two spins, produces a pair of 
spectra each with two peaks whose heights are proportional to the probability 
differences between pairs of states connected by flips of that spin. The factor 
relating the peak heights to the corresponding differences in the probabilities of 
the states can be determined from spectra collected on the pseudo-pure ground 
state, after which it is straightforward to convert the differences into the corre- 
sponding absolute probabilities by linear least squares, subject to the constraint 
that their sum is unity. Wc shall encounter field gradients again in the next sec- 
tion, when we show how they can also be used to implement precisely controlled 
decoherence models. 

Thus the overall experiment consists of collecting ten spectra, as follows: 

1. Prepare the state \P, by first preparing the pseudo-pure ground state | 00) 
using one of the previously described methods, and then transforming it by 
Q(6)P(4>). 

2. Use a selective radio- frequency pulse to apply the rotation R to those spins 
on which measurement B is to be performed. 

3. Use a z-gradient to dephase the transverse components of the resulting den- 
sity operator. 

4. Apply a readout pulse to one of the spins, and collect the corresponding 
spectrum; repeat steps 1-3 and then do the same for the other spin. 

5. Repeat steps 1 - 4 for each of the four combinations of measurements AA, 
AB, BA and BB on the two spins. 

6. Collect two additional amplitude calibration spectra by applying soft readout 
pulses to each spin in the pseudo-pure ground state. 

These experiments were performed on a Bruker 400 MHz spectrometer using the 
two spin i nuclei in 13 C-labeled chloroform. 

The spectra obtained from steps 1 - 5 of this experiment are shown in Fig. 3. 
The probabilities derived from these peak heights, and the residual (square-root 
of the sum of squares of the deviations of the data points from the corresponding 
fit) associated with each, are shown in the table below. It follows that Bell's 
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inequality is violated by 



^AB + ^AA + ^BA_^BB 

= 0.029 + 0.001 - 0.002- 0.097 = -0.069. 

A rigorous error analysis is not possible, because the dominant errors in NMR 
spectra (e.g. RF field inhomogeneity) are not statistical. If we nevertheless take 
the mean RMS residual (half the total residuals shown in the table) of 0.0065 
as an estimate of the errors and assume they are independent between spectra, 
the expected error in the sum of these four numbers is only 0.013, so that this 
violation of Bell's inequality appears significant. 

Nonetheless, as stressed in a recent preprint by Braunstein, Caves, Jozsa, Lin- 
den, Popescu and Schack [11], such experiments on weakly polarized pseudo-pure 
states cannot actually disprove the existence of "hidden variables" associated 
with the spins of the individual molecules. This is because the vast majority of 
a weakly polarized density operator is contained in its identity component, and 
there are many different ensembles of uncorrelated spin states whose net density 
operator is the identity. Hence the noise from the identity component dominates 
the statistics of observations on the ensemble, which are therefore consistent with 
microscopic interpretations in which only uncorrelated states are present with 
nonzero probability. Indeed, if one were to pull the molecules out the pseudo- 
pure sample used in the above experiments one at a time, break them apart, 
and perform the measurements A and B with a Stern-Gerlach apparatus, the 
frequencies of the four combinations of outcomes would all be very close to 1/4, 
and would not violate Bell's inequality. Thus, our apparent violation vanishes 
when the whole ensemble is taken into consideration. 

To see more precisely why the above experiments fail to disprove the existence 
of hidden variables, we first note that a pure state | £ ) ( £ I i s canonically asso- 
ciated with any given pseudo-pure density operator p = (1 — (5)/2 Ar + (5|£)(£|, 
which is distinguished mathematically by the fact that | £ ) is the eigenvector 
corresponding to its sole nondegenerate eigenvalue. We further recall (see Eq. 
(44)) that the traceless part of the pseudo-pure density operator p transforms 
identically to that of the corresponding pure state | £ ) ( £ | under unitary opera- 
tions, while the identity component transforms trivially, and also that p produces 
exactly the same NMR spectrum as would | £ ) ( £ | up to its overall amplitude 
(since the identity component of any density operator does not contribute to 
the signals observed by NMR) . Thus the unitary dynamics of the observables in 
NMR experiments on pseudo-pure states are, for all practical intents and pur- 
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poses, indistinguishable from the same experiments on a (smaller) ensemble in 
the corresponding pure state | £ ) ( £ | . 

It follows that NMR experiments on pseudo-pure states are necessarily con- 
sistent with (though not proof of the reality of) a microscopic interpretation 
of the ensemble in which those molecules contributing to the observations are 
all in the same pure state | £)(£ |, while the remaining (and large majority of 
the) molecules are in completely random states with a net density operator of 
1/2 N . In deriving a violation of Bell's inequality from our measurements above, 
we required that the fractions of molecules in the four diagonal states sum to 
unity, so that they could be identified with the probabilities of those states in an 
unidentified sub cnsemble in the pure state | ip ) ( ip | . Implicitly, therefore, this mi- 
croscopic interpretation of the ensemble was assumed in deriving the violation. 
As explained above, however, the large identity component in the corresponding 
pseudo-pure density operator guarantees that many other ensembles could be 
found with the same net density operator, so that a microscopic interpretation 
in terms of a single well-defined subensemble in the pure state | ip ) ( ip | is not 
physically justified. In fact, the fundamental limits on the amount of information 
that can be extracted on an unknown quantum state even by strong measure- 
ments prevents us from ever knowing if any molecules of our pseudo-pure sample 
exist in or near the corresponding pure state | ip ) ( ip | at all [43] . It is for this rea- 
son that our apparent violation of Bell's inequality fails to disprove the existence 
of hidden variables. 

This ambiguity in the microscopic interpretation of liquid-state NMR exper- 
iments not-with-standing, quantum physics indicates that a psuedo-pure spin 
state, subjected to the same electromagnetic fields as a true pure state, will 
undergo the same unitary transformation. In addition, applying a z-gradient to 
an NMR ensemble renders unobservable the same transverse phase information 
that would be destroyed on performing strong measurements along the z-axis 
on all the spins in the ensemble. Finally, existing experiments relying upon true 
pure states and strong measurements provide direct evidence against hidden 
variable theories (see e.g. Rcf. [2, 10]). Given this background knowledge of the 
underlying physics, our experiments indirectly imply that the pure state | ip)( ip \ 
would violate Bell's inequality. More generally, the ambiguity in the microscopic 
interpretation of liquid-state NMR experiments in no way detracts from their 
utility as a means of studying the dynamics of information contained in either 
pseudo-pure or (by inference) true pure states, even in significantly more com- 
plex spin systems that would be difficult to study by other means. To further 
emphasize this fact, we will now describe NMR experiments we have performed 



31 



which demonstrate quantum error correction using pseudo-pure states. 



6 Quantum error correction by NMR spectroscopy 

The error correcting code we have chosen to illustrate by NMR is well-known in 
the field [35], and uses two ancilla (labeled 2 & 3) to encode the state of a data 
spin (labeled 1). Letting S 211 and S 311 be c-NOT's, and R^ 3 = exp(-«f (I y + 
I y + iy)), the encoding operation proceeds as follows: 

(a|0>+/?|l»|00>/>/2 ^^^l 3 a | +++ ) +/3 |_-_) 
(where|±±±) = (|0)±|l))(|0)±|l))(|0)±|l))/>/8) 

Decoding consists of applying the inverse operations in the reverse order, which 
acts on the states obtained by single sign-flip errors as follows: 

a\++-)+0\ — +) H0) + /?|l))|01)/V2 

a\ +-+ ) + (3\ -+- ) ) + (3\ 1 ))| 10 )/V2 (70) 

a\-++)+f3\+—) R ^° S ^ S ^ Hl)+/3|0))|11)/V2 

It follows that a Toffoli gate T 1 ^ 23 , which flips the data spin conditional on the 
ancillae being in the state 1 11 ), will correct a sign-flip error in the data spin and 
leave it alone otherwise, even if an error occurs in the ancillae. 

In practice, errors in quantum computers are not expected to be single sign- 
flips, but rather small random phase errors which cumulatively result in decoher- 
ence. Nevertheless, we can show that the ability to correct sign-flips implies the 
ability to cancel the effect of such phase errors to first order. Random phase er- 
rors correspond to the propagator exp(— t(x 1 I z + X 2 I 2 +X 3 -^z))' where x 1 , x 2 , x 3 
are random variables, which acts to first order on the encoded state as: 

exp(-*( X 1 / z 1 + x 2 / z 2 + X 3 / z 3 ) ("I +++> + P\ —)) 
w (a\+++)+(3\— ))-i X i(a\-++)+(3\+— )) (71) 
- i X 2 {a\ +-+ ) + 0\ -+-}) - ix 3 H ++- ) + I3\ — +)) 

Since decoding and the error-correcting Toffoli gate are likewise linear, it follows 
that the first-order effects of phase errors are cancelled as claimed. Note this 
argument makes no assumptions concerning the correlations among the errors! 
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Experimental results demonstrating these expectations have recently been 
published [15]. In the following, we shall present a more detailed explanation 
of how the error correction works using the product operator formalism, along 
with selected experimental data illustrating and validating this explanation. We 
shall assume that the data spin is in one of the states 1 (unpolarized), J x , /J or 
l\. Although these are mixed states, each consists of an incoherent sum of pure 
states, e.g. 2l\ = |0)(0| — 1 1>< 1 1, so if error correction works on these pure 
states, by linearity it will also work on the mixtures (and vice versa). In these 
terms, a complete set of initial states p A for error correction are: 



E% E\ 

i\e\e\ 
i)e\e\ 
i\e\e\ 



> = p\E\E\ = p A 



(72) 



The corresponding states p B to which they are mapped by encoding are: 
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+ i 1 i 2 i 3 

1 y z z 


y y y 



(73) 



J(/i+7 x 2 + / x 3 )+/ x 1 J x 2 / x 3 



We note the last three states in Eq. (72) can be prepared (with a 50% loss of 
polarization) from the average of twice Eq. (51) with 



^Diag(3, 1,1, 1,-3, -1,-1,-1) 
<-> 1 L(/ z 1 (3 + 2/ z 2 + 2I 3 Z +4I 2 Z I 3 )) (74) 
= \(E\-E l _)(E\E 3 + + \). 



In liquid-state NMR, decoherence occurs principally through the randomly 
fluctuating external magnetic fields B z along the z-axis at each spin k. The effect 
of these fields is most simply described in the spherical product operator basis 
1, l\ and I± = J x ± ily , as opposed to the Cartesian basis used up to now. 
The products of these basis elements can be shown [19] to decay exponentially 
at rates proportional to the mean-square field (B* ) 2 for J* (as well as I±I*, 
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Fig. 4. Experimental NMR data illustrating the decay of each of the product operators 
l\, H\l\, as functions of the time allowed for gradient diffusion (see 

text), together with the least-squares fits to their logarithms. The single and triple 
quantum coherences in (negative curves) have been plotted and fit separately. 

The sum of these data and of the fits are also shown (topmost curve), which illustrates 
that error correction cancels the decay of the encoded state to first order as expected. 



I k ± I f z r z n ), and to 



for I k + I e _ & 1-1+, 
for I k + f + kl k _l e _, 



and (B* + Bi + B™) 2 for I k + I l + I™ & I k _I l _I m . 



(75) 



These products are referred to as sing le (SQC1: I k ±), zero (ZQC: I k ±I%), dou- 
ble (DQC: I±I±), three-spin single (SQC3: I k ±I e ± I™, etc.) and triple (TQC: 
I±I±I±) quantum coherences, respectively. 

We shall consider two extreme forms of decoherence. In the first, the fields at 
the different spins are uncorrelated, and hence the random variables \ k can be 
assumed to be identically distributed and independent. In the second, they are 
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assumed to be totally correlated. By Eq. (75), the relative rates of decoherence 
in these two cases are: 



ZQC SQC1 SQC3 DQC TQC 
Uncorrelated: 2 13 2 3 

Totally Correlated: 114 9 



(76) 



Decomposing p B into a spherical basis, multiplying by decaying exponentials 
with the above rates normalized by the SQC1 decay rate t, and returning to the 
Cartesian basis gives 



Pc 



\ + {llll + llll + llll)e-^ 
l\l 2 I 3 + Z y Z z 2 Z 3 + I)l*l\)e-^ 
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in the uncorrelated case, and 



(77) 
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in the totally correlated case. The decoding operation converts this to 



Pd 



\ + (ll + ll + llll)e-^ 
{\llll + \llll + llllll)e-^ + ±ll 
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(79) 



in the uncorrected case, and 
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in the totally correlated case. This is clearly getting a little messy, and it gets 
much worse after the Toffoli gate! Therefore, we shall only present the partial 
trace over the ancillae after applying the Toffoli, which is 



Pi 



1 

'y( 



y \2 
3 



-t/T 



l„-3t/r\ 



(81) 



g-t/T _ l e -3t/T) 



3G 



in the uncorrelated case, and 
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(82) 



in the correlated. The slope of these curves at t = is zero in all cases, as 
expected. 

In order to demonstrate these results by NMR solution-state spectroscopy, 
a precise implementation of the above decoherence models is needed. This was 
achieved by combining gradient methods with molecular diffusion. In these meth- 
ods, a magnetic field gradient is created along the z-axis; as previously described, 
this dephases the transverse (xy) magnetization. More precisely, a field gradient 
causes the transverse magnetization to precess at rates which depend linearly 
on its z-coordinate, thereby winding it into a spiral about the z-axis whose av- 
erage transverse magnetization is essentially zero. The gradient is turned off for 
a given time interval t, during which diffusion of the molecules along z blurs the 
spiral. The gradient is then reversed, causing the magnetization to refocus and 
so create an "echo" . Because those molecules which have moved now precess 
at a different rate, their magnetization is not rcfocussed, so the magnitude of 
the echo decays exponentially with t. Because all the spins in each molecule are 
subject to the same change in field, this constitutes a true implementation of 
the totally correlated model. By using rcfocusing 7r-pulscs between gradients, it 
is also possible to dephase each spin separately, thereby implementing the un- 
correlated model. At this time, however, we have collected and processed data 
only for the p\ — l\ state with the totally correlated model. 

Although it is possible to prepare the state I Z E + E + as noted above, we have 
chosen to illustrate the above analysis by preparing the states J z , 2J Z 7 Z , 2J Z J z 
and in four separate experiments, each using sixteen different decoher- 

ence times t. Because the SQC and TQC contributing to 4J Z refocussed at 
different times, this further enabled us to follow their evolutions separately. The 
results of these experiments are plotted against the time t in Figure 4, along 
with the corresponding logarithmic fits. It may be seen that the sum of the data 
and of the fits thereto (also shown) do indeed exhibit a near-zero initial slope, 
in accord with the above calculations. Our published report [15] includes the re- 
sults of further experiments (performed by E. Knill and R. Laflamme) with the 
natural and far more complicated decoherence processes that occur in solution. 
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These arc more difficult to interpret, but are nevertheless consistent with the 
state preservation expected from error correction. Additional experiments and 
more detailed calculations are in progress. 

While a method of inhibiting decoherence (T2 relaxation) during NMR pulse 
sequences would be highly desirable, there are strong reasons to doubt that 
quantum error correction will be useful in this regard. First, the ancillae must 
be placed in a pseudo-pure state, which as we have shown above entails a loss 
of 50% of the signal for each "data" spin; this is more than is recovered by error 
correction. In addition, the ancillae must be returned to a pseudo-pure state 
uncorrelated with the state of the data spin(s), or else "fresh" ancillae in such a 
state must be continuously available, in order to inhibit decoherence over an ap- 
preciable period of time by the repeated correction of errors. Nevertheless, we feel 
that the basic idea underlying error correction of preparing multiple quantum 
coherences, allowing them to decohere, and then mixing them so as to determine 
their relative rates of relaxation, may be of considerable use in NMR studies of 
the statistics of molecular motion. This in turn is one of the most important ap- 
plications of NMR spectroscopy. Conversely, whereas NMR spectroscopists have 
previously used their methods solely to unravel the secrets of naturally occurring 
systems, it now appears possible to use these same methods to engineer artifi- 
cial systems in which the basic principles of quantum information processing, in 
particular the emergence of the classical world through decoherence [23] , can be 
studied in unprecedented detail. 
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